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FOREWORD 


This  report  was  prepared  by  the  Polytechnic  Institute  of  Brooklyn 
under  two  contracts.  The  experimental  work  and  data  reduction  were 
performed  under  USAF  Contract  AF  33( 6l 6)-5944;  the  theoretical  rnalysis, 
interpretation  of  the  data)  and  preparation  of  the  report  were  performed 
under  USAF  Contract  AF  33{6l6)-766l.  The  contracts  were  initiated  under 
Project  No,  7364,  "Experimental  Techniques  for  Materials  Research," 
Task  No.  73652,  "Intense  Thermal  Energy  Transfer  Into  Materials,"  and 
Project  No.  7064,  "Research  on  Aerodynamic  Flow  Fields,"  Task  No. 
70169,  "Thermal -Aerodynamic  Characteristics  at  Hypersonic  Mach 
Numbers,"  respectively.  The  work  was  administered  by  the  Materials 
Central,  Deputy  of  Advanced  Systems  Technology,  Aeronautical  Systems 
Division  and  the  Aeronautical  Research  Laboratory,  Air  Force  Office  of 
Aerospace  Research,  with  Mr.  Hyman  Marcus  and  Lt.  John  D.  Anderson 
acting  as  project  engineers,  respectively. 

The  report  covers  work  conducted  from  February  I9  60  to 
February  1961.- 

*  This  report  was  prepared  by  Dr.  Robert  J.  Cresci,  Research 
Assistant  Professor  of  Aeronautical  Engineering  and  Dr.  Paul  A.  Libby, 
Professor  of  Aeronautical  Engineering  and  Assistant  Director  of  the 
Aerodynamics  Laboratory. 

The  authors  are  pleased  to  acknowledge  the  valuable  assistance 
of  the  hypersonic  and  computing  staffs  of  PIBAL  in  carrying  out  this 
research,  and  of  Professor  Lu  Ting  for  suggestions  concerning  the  analysis 
of  the  shock  stand-off  distance. 


ABSTRACT 


The  influence  of  localized  mass  transfer  at  the  nose  of  a  slender 
con-3  under  hypersonic  flow  conditions  has  been  studied  by  experimental 
and  theoretical  means.  Two  gaseous  coolants,  nitrogen  ang  helium  are 
injected  through  a  porous  plug  subtending  a  half  angle  of  30  .  The  effect 
of  the  mass  transfer  on  the  shock  shape,  pressure  distribution,  heat 
transfer  and  transition  are  investigated.  The  experimental  work  involved 
tests  in  the  Mach  number  8.  0  tunnel  at  PIBAL,  The  theoretical  analysis 
involved  a  study  of  the  effect  of  mass  transfer  on  the  shock  stand-off 
distance  and  leads  to  an  inviscid  flow  parameter  permitting  the  experi¬ 
mentally  determined  shock  shape  and  pressure  distribution  to  be  extra¬ 
polated  to  other  than  test  conditions  and  to  other  coolant  gases.  There 
is  obtained  the  maximum  value  of  this  parameter  resulting  in  no  signi¬ 
ficant  alteration  of  the  pressure  distribution  on  the  cone  and  thus 
defining  the  flows  in  which  boundary  layer  type  of  similarity  applies. 

Significant  reductions  in  heat  transfer  are  obtained  with  injection. 
Indeed  with  small  amounts  of  helium  injection  the  peak  heating  is  found  to 
occur  downstream  on  the  cone  and  to  be" an  order  of  magnitude  less  than 
would  occur  at  the  stagnation  point.without  mass  transfer.  With  nitrogen 
early  transition  is  found  to  occur  so  that  local  heating  rates  are  actually 
increased  over  those  prevailing  at  the  same  Reynolds  number  without 
injection. 
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I  -  INTRODUCTION 


The  use  of  mass  transfer  for  the  reduction  of  heat  transfer  has  been 
the  subject  of  considerable  •research  in  connection  with  high  performance  . 
propulsion  systems  and  with  hypersonic  flight.  As  a  residt  a  variety  of  heat 
protection  systems,  some  involving  mass  transfer  have  been  developed. 

The  research  reported  herein  pertains  to  the  use  of  localized  mass 
transfer  for  the  reduction  of  Keat  transfer  in  the  nose  region  and  on  the 
downstream  conical ’portion  of  a  slender ’cone.  “The  syst’em.  considered  is 
shown  schematically,  in  Figure  la.  A  relatively  slender  cone  has  a  spherical 
nos'e  made  of  a  porous  material  through  which  a  gaseous  coolant  is  injected. 
The  coolant  has  two  effects  on  the  heat  transfer;  it  reduces  and  absorbs  the 
heat  transferred  to  the  spherical  region  and  reduces  the  heat  transferred  to 
the  downstream  conical  portion  of  the  body.  In  addition,  "the  coolant  can 
alter  the  shock  shape  and  the  pressure  distribution  o"n  the  body. 

This  localized  mass  transfer  system  is  conceptually  similar  to  the 
system  studied  in  reference  1.  There  a  hemisphere  cylinder  with  a  porous 
plug  subtending  an  angle  of  roughly  10  with  the  axis  of  symmetry*  was 
considered.  The  gaseous  coolant  injected  through  the  plug  reduced  the 
heat  transfer  at  the  stagnation  point  and  on  the.  downstream  portion  of  the 
hemisphere.  The  coolants  tested  were  helium,  nitrogen,  argon  and 
krypton. 

In  considering  the  similarity  conditions  for* these  localized  mass 
transfer  systems  the  present  authors  in  reference  1  distinguish  between 
two  types  of  flows;  those  in  which  the  coolant  is  introduced  into  a  classical 
boundary  layer  without  significant  alteration  of  the  inviscid  flow  as  shown 
in  Figure’ lb,  and  those  in  which  the  entire  flow  field  is  influenced  by  the* 
coolant  flow  as  shown  in  Figure  lc.  For  the  first  type  of  flow,  similarity 
conditions  based  on’ the  boundary  layer  equations  are  applicable  and  can 
be  employed  to  extrapolate  tes.t  data  to  other  Reynolds  numbers  and  to 
other  Mach  numbers  if  the  pressure  distribution  is  Mach  number  insensi¬ 
tive.  On  the  other  hand,  flows  of  the  second  type  lead  to  stringent  simili¬ 
tude  requirements;  indeed,  no  extrapolation  of  test  data  to  other  Reynolds, 
and  Mach  numbers  is,  in  general,  possible.  These  flows  involve  viscous 
layers  which  are  relatively  thick  compared  to  the  shock  layer  and  to  the 
cylindrical  radius.  Consequently,  the  shock  shape  and  pressure  distri¬ 
bution  change  as- either  the  rate  of  mass  transfer,  or  the  Reynolds  number, 
or  the  Mach  number  is  changed. 

•The  experimental  results  in  reference  1  indicated  that  for  the 
porous  material,  flow  conditions,  and  rates  of  coolant  mass  flow  employed 
therein,  the  boundary  layer  was  laminar  and  the  inviscid  flow  field  was 
essentially  unaltered  at  even  the  highest  rates  of  mass  transfer;  therefore 
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POROUS  REGION 


the  flow  was  of  the  first  type.  Accordingly,  the  similarity  parameters 
developed  in  reference  1  could  be  used  to  extrapolate  the  test  results  either  to 
flight  or  to  other  test  conditions  provided,  of  course,  that  the  flow  remains 
laminar. 

*  •  • 

•  • 

.  For'the  localized  mass,  t  ran  sfe't  system  disdussed  herein  the  same 

considerations  with  respect  to  flow  similarity  will  apply.  For  s'ufiiciently 
low  rates  of  coolant  flow,  the  inviscid  flow  field  will  be  unaltered  by  the 
injection  and  similarity  parameters  prtesented  in  reference  1  are  applicable 
as  follows:  Consider  two  geometrically  similar  bodies  as  shown  schem¬ 
atically  in  Figure  la,  with  identical  inviscid 'flows  in  terms  of  the  distri¬ 
bution  of  p  =  p(x),  u'=  u(x),  etc.  Assume  that  the.chemical  composition's 
of  the  coolant  and  of  the  external  gas.  stream  a^e  fixed  and  equal  in  the 
flows.  Assume  further  that  the  enthalpy  ratio  h'c  (coolant-tp-free- stream 

stagnation)  and  the  distribution  of  temperature  ratio  0^  (wall-to-free- 

stream  stagnation)  are  the  same  in  the  two  flows.  Then  the  rates  of 
laminar  heat  transfer  to  the  impermeable  surface  on  the  two  bodies  are 
related  by  the  similarity  parameters  ***. 


N  ='N  (N  , he 6  ,  k) 


(1) 


I  where  N  and  N  are  respectively  mass  and  heat  transfer  parameters 
l  a 

•  defined  as  follows  ;*  *  ... 


N  .  mc/ v„  NR! 


(2)- 


and 


N  s  N  *  /Nj 
a  Nu  •  R 


=  CqR  (c  )  7(h  -h>  '  V(p'hjR./us)(ps  /p-h  )lyl  (3) 


o  p.s  s  w  8  -  rs  s.  o'  s  “s  rs  s 
r  e  e  e  •  e  e  eeee 


•  • 

This  result  clearly  permits  heat  transfer  data  obtained  at  one  Reynolds 

number  (N^)  to  be  extrapdlat'ed  to  other  flow  conditions  provided  N  ,  hc> 


6  and  x.  are  unaltered, 
w  1 


*  Consider  the  theoretical  determination  of  the  boundary  layer  char¬ 
acteristics  of  flows  covered  by  these  similarity  requirements.  If  such  a 
determination  were  possible,  the  functional  relation  expressed  by  Eq.  (1) 
could  be  evaluated  by  theoretical  analysis.  Two  regions  along  the  surface 
can  be  identified;  for  x<x.  the  boundary  layer  involves  mass  transfer, 

•  while  for  x>x,,  the  surface1  is  impermeable.  Now  it  was  found  in  refer- 
‘ence  1  that  the  rates  of  injection  required  to  obtain  a  significant  down¬ 
stream  effect  from  stagnation  point  mass  transfer  were  considerably 
.higher  than  those  required  to  reduqe  substantially  the  heat  transfer  to  the 


•  *  ( 

porous  surface  itself.  Indeed,  the  values  of  mass  transfer  per  unit  ^rea 
per  iinit  time  (pv)  employed ‘in  those  experiments  were  found  to  be  up  to 
20  or  30  times  thc$3e  usually  considered  in  theoretical  analyses  of  trans¬ 
piration  cooling.  It  would  be  expected  a  priori  that  similarly  high  rates 
of  mass  transfer  would  be  required  in  th'e  configuration  considered  here 
for  essential  reductions  in  the  heat  transfer'to  the  downstream  conical 
section.  ' 


In  reference  2  a  theoretical  analysis  bf  the  stagnation  point 
boundary  layer  with  large  rates  of  air-to-air  injection  was  carried  out. 

It  was  found  that  the  boundary' layers  could  be  characterized  by  an  iso¬ 
thermal  shear  layer  close  to  the  surface  and  a  relatively  thin* mixing  zone 

•  adjusting  this  inner  flow' to  the  'external  flow*  values.  This  result  <  iggested 
an  integral  method  which  was  shown  to  lead  to  accurate  values  for  the 
gross  boundary  layer  properties.  *  This  .characterization  of  the  stagnation 

•  point  boundary  layer  also  suggests  the  idealization  which  will  be  employed 
herein  for  the  study  of  the  alteration  of  shock  stand-off  distance' with 
injection;  accordingly,  the  boundary  .layer  is  .replaced  by  an  isothermal 

shear  layer  which  is  matched  to  the  external  flow  by  a  slip  surface. 

• 

It  appears  then  that  the  analysis  of  the  first  region  (0  <  x  <  x.)  for 
large  rates  of  injection  can  be  carried  out  according  to  either- of  two 
approximate  methods.  The  stagnation  point  analysis*for  large  rates  of 
.  .  injection 'can  be'considered. to. apply  approximately  over  the  entire  range 
•  of  x.  In  general,  the  requirements  of  boundary  layer  similarity,  especially 
‘with  respect  to  the  distribution  of  the  rate  of  mass  transfer.,  will  not  be 
exactly  satisfied;  however,  a  stagnation  point  analysis  should  be_sufficiently 
•accurate  fdr  most  purposes.  A  second  method  of  analysis" would  be  to 

.  apply  an  integral  method  as  in  reference  2  to  the  development  of  the 
.  boundary  layer  over  the  entire  range  of  x  starting  with  stagnation  point 

•  •”  values.  *  * 

•  •  • 

•  • 

The  analysis  of  the  boundary  layer  in  the  downstream  impermeable 
region  (x  :>  x.)  requires  the  treatment  of  the  initial  value  problem  of 
.  boundary  layer  theory.  At  x  =  x.  the  velocity,  stagnation  enthalpy  and 
concentration  profile's  are  given xby  the  analysis  of  the  porous  region; 
the  downstream  development  of  the  boundary  layer  is  desired.  This 
problem  can  be  treated  by  integral  methods  as  in  references  3  and*4, 
by  finite  difference  methods  as  in  reference  5,  or  by  the  method  of 
"strips",  recently  published  by  Pallone  (reference  6)  and  more  recently 
described. by  Dorodnitsyn.  ** 


35C  •  • 

•  The  analysis  of  reference  2  is  at  present  being  extended  to  the  case  of 
*  •  *^elium-to-air  injection,'.' 

At  the  Second  International  Conference  of  the  Aeronautical  Sciences, 
Zurich',  Switzerland,  September  11-  16,  1960. 
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In  this  report  there  are  presented  the  results  qf  an  experimental # 
investigation* of  a  configuration  close  to. that  shown  schematically  in 
Figure  la.  The  tests  were  carried  out  at  a  Mach  number  of  8.  0  in  the 
hypersonic  facility  of  PIBAL  (cf.  references  7  and  8).  Two  coolants,  . 
helium  and  nitrogen,  were  injected  over  sufficiently  wide  ranges  of  mass 
rates  and  Reynolds  numbers,  such  that  both  types  of  flow  described 
herein  prevailed. 

• 

Three  types  of  information  were  obtained  from  these  tests;  the 
first  pertains  to  the  alteration  of  the  shock  shape  and  pressure  distribution 
by  the  injection,  A  theoretical  analysis  of  the  influence  of  mass  transfer* 
on  the  shock  stand-off  distance  is  carried  out  and  shown  to  yield  results" 
in  good  agreement  with  experiment.  The  well-known  approximations  of  . 
a  spherical  shock  with  a  constant  density  shock*layer  are  employed  al*ong 
with  an  idealization  of  an  inner  layer  formed  by  the- injected  fluid;  namely,  . 

*ah  isothermal  shear  layer ‘with  vorticity.  .  This  theoretical  model,  which 
inapplicable  to  both'nitrogen  and  Helium  injection,  suggests  a  combination 
of  N^*  andNp,^  i.  e.  ,  a  new  inviscid  flow  parameter  Ng.  From  the 

available  pressure  data,  the  maximum  value  of  N  for"  no  alteration  in 
pressure"  distribution  is  obtained.  *  •  3  •  •  • 

■  .  .•••**. 

.  .  The  second  type  of  data  presented  in  the  report  is  the  heat  transfer 
to  the  downstream  conical  surface*.  Laminar,  transitional  and  turbulent 
heat  transfer  rates  are  obtained  and  compared  to  the  available  theories 
for  zero  mass  transfer.  “Significant  reductions  in  .laminar  heat  transfer 

.result  from  thie  mass..trhnsfer  but  earlier  transition  is  found  to  occur. 

•  •  •  %  • 

.  *  .  • 

.  Finally,  the  third  type  of  data  presents  the  alteration  of  transition 
behavior  with  Reynolds  number  and  maqs  transfer.  Tiiese  data  are 
limited;  additional  tests  would  be  required  for  reasonable  knowledge  of  * 

.  transition  behavior. 

•  • 


H.  MODEL  DESIGN,’  TEST  PROCEDURES  AND  TEST  CONDITIONS  •  . 

The  modeloused  in  the  tests  is.  a  blunt-nosed  cone  with  a  semi¬ 
vertex  angle  of  20  .  The  injection  region,  consisting  of  a  spherically 
shaped  porous  cap,  exten’ds  over  a  half  angle  of  30*.  The  details  of 
the  model  design  are  shown  in  Figure  2  along  with  the  the  thermocouple  and 
pressure" tap  locations.  Figure  3  presents  photographs  of  the  model  on 
the  support" prior  to  installation  in  the  wind  tunnel. 
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FIG.  (2)  DETAILS  OF. MODEL  DESIGN  AND  INSTRUMENTATION 


,  •  • 

The  wind  tunnel  tests  presented  herein  were  performed  at  g  Mach 
number  *of  8.  0,  at  a  stagnation  temperature  of  approximately  1700.  R, 
and  at  free  stream  stagnation  pressures  varying  between  150  and  600  psia, 
Schlieren  pictures  were  taken  during  each  test.  The  maximum  mass  flow 
injected  at  the  model  nose  was  0.  01  lb/sec.  ‘for  helium  and  0.  03  lb/sec. 
for  nitrogen.  In  terms.of  the  nondimens ional  similarity  parameters 
defined'in  the  Introduction,  the  Reynolds  number  variatipn  corresponds 
to  0.  7  <  N„  x  10_4,<  2. ’5,  and  the  mass  transfer  parameter  variation 
N  <‘  7?" for  helium  and  N  *<213  for  niFro'gen".  The  coolant  temperature 
i  —  •  •  l  —  *  .  o 

.  remained  essentially  constant  for  all  tests  and  was  approximately  530  R. 

On  the  impermeable  surface  downstream  of‘the  injection  region,  the 

average  wall  temperature  was  560  R  for  all  th’e  tests;  thus,  the  heat 

transfer  data  correspond  to  an  isothermal  surface  with  0  -  0.  33.  ‘The 

model  was  initially  polished  to  a  surface  roughness  of  less  than  4  micro- 

inches.  As  the  test  series  proceeded,  surface  erosion  occurred;  however, 

the  model  surface  was  repolished  so  that  the  maximum  roughness  is 

believed  to  be  in  the  neighborhood  of*30  to  40jjl  for'all  tests. 

•  •  • 

•  ^  § 

Pressure  data  were  obtained  by  differential  transducers  and  re¬ 
corded  on  single  channel  recording  potentiometers.  In  order  to  minimize 
the  response  time,  the  reference  chamber  of  the  transducers  was  evacuated 
to  a  pressure  slightly  below  the  static  pressure  in  the  tunnel  during  the  run. 
W.ith  this  technique  the"  steady  state  indication  of  *pressure*was  achieved  in 
about  5  seconds,  permitting  relatively  short  running  times. 


* 

The  heat  transfer  data  were  obtained  by  a  transient  technique 
utilizing,  a  thin  skinned  model  at  an  essentially  uniform  initial  temperature 
The  local  heat  transfer  rate  can  be  related  to  the  initial  timewise  deriva- 
.  tive  of  the  surface  temperature;  the  details  of  the  analysis  and  the  method 
of  thermocouple  installation  are  discussed  in  detail  in  references  1  and  9, 
fpr  example.  *  .  *  . 


The  mass  transfer  rate  through  the  porous  plug  was  determined  . 
by  two  independent  measurements.  A  flow  meter  in  the  coolant  supply 
line  provided  a  direct  measurement  of  the  coolant  flow;  the  measured 
pressure  drop  across  the  plug  also  determines  the  mass  flow  rate  according 
to  the  following  relation: 


0,  55  r  2 

•m  =  CW.  LP 
c  c  r  c 


-  V3 


anO.55 


(4) 


The  constant,  C,  was  obtained  by  calibrating  the  porous,  plug  over  the 
•entire  mass  flow  range  prior  to  installation  in  the  tunnel.  The  mass, 
transfer  rates  as  determined  by  the  two  measurements  differed  By  a* 
maximum  of  10%  in  all  tests.  *  Co’olant  temperature  was  measured  by  an 


i's  of  interest  to  note  that  lor  th'e  particular  porous  material  used  in 
these  tests  p  3  >>p  3  so  that  m  is  essentially  constant  with  x.  Thus 
the  agreement  between  the  calibrations  with  zero  and  hypersonic  free 
stream  velocity  is  not  surprising.  •  . 


open  tip  thermocouple  located  along  the  coolant  chamber  canterline.  . 

• 

•  The  porous  cap  was  formed  by  coining  a  flat  sheet  of  porous  stain¬ 
less  steel.  Initially,  the  mean  pore  opening  was  65  microns;  however, 
as  a  result  of  the  forming  process  the  openings  on  the  outer  surface 
increased  slightly.  To  as.certain  qualitatively  the  effect  .this  would  have 
on  the  distribution  of  inje'ction  velocity,  the  spray  pattern  with  water  as 
a  coolant  was  observed.  Various  plugs. were  tested  until  one  was  obtained 
with  a  relatively  uniform  spray  distribution.. 

•• 

The  accuracy  of  the  measurements  is  estimated  to  be  as  follows: 

the  overall  error  in  the  pressures  on  the  cone  surface  is  ±  .  02  psi.  The 

mass  flow  pf  the  coolant  gas  is  in  error  by  a  maximum  of  ±  5%  over  the 

entire  flow  range.  The  heat  transfer  results  involve  many  parameters 

affecting  the  final. dataj  an  o.veralh  estimate  of  the  error  is  believed  to  be 

±  *10%  in  terms  of  the  Nusselt  number.  . 

*  *  *  . 


HI.  PRESENTATION  AND  DISCUSSION  OF  RESULTS 

*  .  • 

•  « 

•  •  • 

•  . 

The  data  obtained  in  the  tests  present  three  different  effects 
of  the  mass  transfer;  on  the  pressure  distribution  and  shock  shape,  on 
the  heat  transfer  to  the  cone ,  and  on  transition.  The  data  pertaining 
thereto  are  discussed  iri  the  following  sections. 

Pressure  Distribution  and  Shock  Shape 
- . - - - -  • 

.  Figure's  4a-4c  and  5a-5d  present,  respectively,  the  pressure 
data  and  typical  schiieren  photographs.  In  Figure  4a  the  pressure  dis- 
tribution.for  zero  mass  flow  is  presented  as  a  line  representing  the  mean 
of  the  data  points.  Also  shown  is  .the  theoretically  predicted  pressure 

distribution  fbr  M  =8.  0  with  T  =5000  R  obtained  from  reference  10,  and 

co  .  s 

•  •  • 

the  sharp  cone  value  for  a  calorically  perfect  gas  from  reference  11,  for 
example.  The  data  points  shown  in  Figure  4a  along  with  the  values  of 

N  and  N_  pertaining  thereto  correspond ’to  no  alteration  of  the  pressure 
.  1  R 

distribution  within  experimental  accuracy.  In  Figures  4b  and  4c  there 
are  shown  the  pressure  distributions  which  result  from  relatively  large 
rates  of  nitrogen  and.helium  injection,  respectively.  Again  the  values 
of  N  and  N_  for  each  distribution  are  giverf. 

i  R  , 

From  these  figures  and  additional  schiieren  photographs  the 
following  observations  may  be  made: 

• 

1)  The  experimentally  determined  pressure  distribution  over  the 
entire  cone,  for  zero  mass  transfer  is  within  10%  of  that  theoretically 
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FIG. (4)  PRESSURE  DISTRIBUTION  (a)  NO  ALTERATION 

DUE  TO  INJECTION 


FIG. (4)  PRESSURE  DISTRIBUTION  (c)  HELIUM  INJECTION 


FIG. (5)  TYPICAL  SCHLIEREN  PHOTOGRAPHS  WITH  VARIOUS  MASS.  TRANSFER  RATES 


predicted  in  reference  10,  This  result  is  in  agreement  with  the  generally 
accepted  insensitivity  of  the  pressure  distribution  to  real  gas  effects  for 
stagnation  temperatures  below  roughly  5000  R, 

2)  The  effect  of  large  rates  of  injection  is  to  lower  the  static 

pressure  on  the  surface  of  the  cone. 

•  » 

3)  At  a  value  of  =2.  2  x  104  ,  alteration  of  the  pressure  distribu¬ 
tion  on  the  cone  surface  is  observed  for  N  >  73.  7  for  nitrogen  and 

N  >  24.  9  for  helium.  1 

l  — 

4)  Examination  of  schlieren  photographs  with  increasing  injection 
rate  and  decreasing  Reynolds  number,  indicates  that  the  shock  stand-off 
distance  is  increased  and  that  this  alteration  occurs  before  an  observable 
change  in  pressure  distribution  occurs. 

•  '  In  order  to  extrapolate  the  heat  transfer  data  presented  here  to 
different  flow  conditions,  the  combinations  of  injection  rate  and  Reynolds 
number  for  which  the  pressure  distribution  is  essentially  unaltered  must 
be  estimated.  Since  the  change  in  pressure  distribution  is  associated  with 
a  radical  change  in  the  shock  stand-off  distance,  a  theoretical  analysis 
leading  to  the  prediction  of  the  dependence  of  shock  stand-off  distance  on 
mass  transfer  has  been  carried  out  and  presented  in  Appendix  I,  This 
analysis  predicts  that  the  shock  stand-off  distance  depends  on  an  inviscid 

l  •  l 

parameter  N  s  m  (0  W  )s  /  [Roa(WcPg  pg  )a  ]  which  can  be  formally 
.  l  6  e 

identified  as  N  s  N  (0  W  /N„W  )S.  *  Accordingly,  there  are  presented 
a  i  w  6  Jtv  c 

•  in  Figure  6  the  variations  with  this  parameter  of  A/R  as  obtained  from 
the  schlieren  photographs.  ^  It  will  be  observed  that  the  correlation  of  the 
experimental  data  for  both  helium  and  nitrogen  is  good,  and  that  i/R  is 
increased*by  an  order  of  magnitude  by  .injection.  Also  presented  in  Figure 
6  are  the  theoretical  predictions  of  A/R  given  in  Appendix  I  and  based 
either  on  a  consistent  theory  or  on  an  empirical  correction  of  the  velocity 
gradient.  It  will  be  noted  that  both  theoretical  predictions  agree  well  with 
the  data  although  the  empirical  correction  l'eads  to  somewhat  better  agree¬ 
ment. 

Since  the  shock  stand-off  distance  was  found  to  be  dependent  on 
N  ,  it  was  of  interest  to'see.  if  the  shock  shape  and  the  pressure  distribu- 
tiSn  on  the  cone  can  be  correlated  in  terms  of  N  .  Comparison  of  Figures 
5c  and  5d  and  of  other  schlieren  photographs  whtrein  N  is  roughly  constant 
although  Nx  and  the  coolant  fluids  are  different,  indicates  that  the  corre¬ 
lation  between  shock  shapes  is  good.  *  • 

•  % 

*’j{C  '  "r  '  '  "  ~  i  ""  ~  " 1  -1-  -  -  'I  ■  ■  "  T  l -  -  -ni---  -  r  —  „  „  , 

Note  that  p.  drops  out  of  N  so  that  it  is  independent  of  viscosity.and 
is  thus  an  inviscid  flow  parameter. 
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FIG  .(6)  VARIATION  OF  SHOCK  STAND-OFF  DISTANCE  WITH  INJECTION  RATE 


,  A  careful  study  of  th'e  corresponding  pressure  distributions  indi¬ 
cated  that  the  experimental  accuracy  was  not  sufficient  to  establish  N  as 
the  correlation  parameter  for  the  pressure  distribution.  The  correlation 
of  shock  shape  with  N  implies,  however,  that  more  accurate  pressure 
measurements  would  confirm  the  correlation  of  pressure  distribution  in 

terms  of  N  . 

3 

Now  consider  the  value  of  N  which  leads  to  an  observable  altera¬ 
tion  of  the  pressure  distribution.  ^ec^use  of  experimental  error  it  is 
difficult  to  select  a  precise  value  therefor;  however,  the  aforementioned 
careful  study  of  the  pressure  distributions  leads  to  the  result  which  is 
shown  in  Figure  7,  and  which  defines  roughly  the  critical  values  of  Ng, 
namely,  l/4<N  <l/3.  Thus,  for  values  of  N  less  than  these,  the  pressure 
distribution  is  unaltered,  the  flow  ia  of  the  boundary  layer  type,  and  the 
parameters  N  and  can  be  employed  in  determining  the  heat  transfer 
for  other  thanHest  conditions. 


Heat  Transfer 

The  heat  transfer  data  for  all  tests  are  presented  in  Figures  8a-8m 
in  terms  of  the  heat  transfer  parameter  N  versus  N  at  each  thermocouple 
location.  Data  for  each  gas  and  the  value  ^>f  N„  for  £ach  test  are  presented 
so  that  the  Nusselt  number*at  each  measuring  station  and  in  each  test  can 
be  explicitly  determined  if  desired. 

.  .  . 

The  interpretation  of  these  heat  transfer  results  is  Complicated  by 
*two  effects;  ’the  alteration  of  the  pressure  distribution  by  injection  as 
related  to  the  parameter  N  and  the  occurrence  of  transition.  If  neither 
’  of  these  two  effects  were  present,  it  would  be  expected  tha.t  N  w;ould 
monotonicaliy. decrease  with  increasing  Nx  at  all  thermocouple  locations 
except  those  far  downstream  where  N_  would  be ’expected  to  be  independent 
of  N.  .• 

1  • 

•  ^ 

.In  discussing  these  effects  it  is  convenient  to  consider  a  fixed 

and  a  fixed. thermocouple  location.  The  observed  reduction  in  static  pressure 
on  the  cone  might  a  priori  be  expected  to  lead  to  further  reductions  in  heat 
.transfer  with  increasing  Nx  .  However,  the  response  of  the  relatively  thick, 
inner  layer  of  injected  gas  to  simultaneous  viscous  shear  and  to  pressure 
gradient  is  complex;  indeed,  «it  will  be  seen  that  for  the  case  of  helium 
injection  involving  a  different  ratio  of  specific  heats,  the  alteration  of 
pressure  distribution  is  accompanied  by  an  increase  in  heat  transfer.  If 
transition  is  associated  with  increasing  Nx,  the  heat  transfer  would  be 
expected  to  increase.  To  assist  in  interpretation,  the  data  in  Figures  8a-8m 


FIG. (7)  CRITICAL  REYNOLDS  NUMBER  VS  MASS  TRANSFER 
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FIG.  (8)  VARIATION  OF  HEAT  TRANSFER  WITH  MASS  TRANSFER. (d)x*4.84 
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FIG.  (8)  VARIATION  OF  HEAT  TRANSFER  WITH  MASS  TRANSFERS i  )*=9.84 
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FIG.  (8)  VARIATION  OF  HEAT  TRANSFER  WITH  MASS  TRANSFERS  I)  5c *12.84 
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have  been  identified  with  a  flag  if  N  >l/3,  i.  e. ,  the  value  of  Ns  resulting 
in  an  observable  change  in  pressuA  distribution. 

Consideration  of  Figures  8a- 8m  indicates  the  following: 

•  1)  For  a  given  Na  the  reduction  in  Na  is  considerably  greater  for 

helium  than  for  nitrogen. 

2)  For_helium,  N  decreases  sharply  with  increasing  N  for 
N1<25  and  for  s<8.85.  l?or  N>25  and  s>8.85,  N  is  relatively 
insensitive  to  N1 .  For  a  fixed  and  increasing  Nx ,  the  heat  transfer 
increases  when  N3iSl/3;  further  increases  in  lead  to  a  second  decline 
in  heat  transfer.  From  a  consideration  of  the  data  for  fixed  Nx  but  varying 
ftp ,  it  is  not  possible  to  determine  if  transition  is  occurring  in  these  tests. 

•  •  • 

3)  For  nitrogen  there  appears  to  be  no  increase  in  heat  transfer 
associated  with  the  alteration  in  pressure  distribution,  i.  e.  ,  with 
Na>l/3.  JHowever,  the  difference  in  the  values  of  Na  for  a  given  Nx  but 
different  implies  that  transition  is  present.  *  Indeed,  there  appears  to 
be  only  one  nitrogen  test  in  which  both  laminar  flow  and  unaltered  pressure 
distribution  occurred.  Transition  is  seen  to  lead  to  heat  transfer  rates 
which  are  higher  with  mass  transfer  than  without. 

•  * 

4}  With  values  of  Nj-  5  for  helium  the  heat  transfer  over  the 
entire  cone  corresponds  to  Na  -  0.08  or  less.  This  should  be  compared 
to  Na=0.84  which  exists  at  a  stagnation  point  without  mass  transfer 
(cf.  references  1,  2  and  9).  Note  that  on  the  porous  regions  the  heat 
transfer  is  zero  fcjr  relatively  Small  values  of  Nx  as  may  be  inferred 
from  reference  2.  Thus  as  a  heat  protection  system  the  configuration 
studied  here  is  quite  effective’.** 

*  *  *  • 

It  is  of  interest  to  consider  two  cross-plots  of  the  data  of ^Figures 
8a-8m.  In  Figure  9  there  are  given  the  •distributions  of  Na  with  x  for 
those  tests  in  which  the  inviscid  flow  was  considered  unaltered  by  injection 
(N^.l/3)  and  in  which  laminar  flow  prevailed.  Also  shown  thereon  is 
the  theoretically  predicted  distribution  of  Ng  for  zero  injection  obtained 
from  a  combination  of  the  analyses  of  references  12  and  13.  The  distri¬ 
bution  of  heat  transfer  relativ.e  to  the  stagnation  point  value  is  computed 
from  reference  12  with  a  constant  isentropic  exponent  of  1.32  and  the 


The  transition  data  will  be  discussed  in  detail  below. 

It  may  be  of  interest  to  note  that  if  R  =0.5  ft. ,  h=50,  000  ft.  and 
V  -  8000  fps,  m  =0.19  lbs/sec  for  I?  =5.0. 

CO  c  *1 
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measured  pressure  distribution}  the  stagnation  point  value  is  obtained 
from  reference  13.  The  stagnation  pressure  behind  a  normal  shock,  was 
assumed  to  prevail  over  the  entire  body:  It  will  be  noted  that  the  agreement 
between  theory  and  experiment  for  zero  injection  is  within  20%.  It  will  also 
be  seen  from  Figure  9-  that  substantial  reductions  in  laminar  heat  transfer 
over  extended  portions  of  the  cone  occur  with  helium  injection.  The 
influence  of  injection  of  helium  is  seen  to  persist  to  some  extent  even 
around  the  junction  of  the  cone- cylinder. 

A  second  cross-plot  of  interest  is  shown  in  Figure  10;  there  the 
distributions  of  Nusselt  number  Nj^  with  x  are  given  for  nitrogen  as  a 

coolant  and  for  a  fixed  Reynolds  number,  =  2.  2(104).  In  addition  to 

the  experimental  data. ,  there  are  shown  the  theoretical  predictions  for 
laminar  flow  as  discussed  above  in  connection  with  Figure  9  and  for  fully 
developed  turbulent  flow.  The  flat  plate- reference- enthalpy  (F.  P.  R.  E. ) 
method  (cf.  reference  14)  was  employed  for *turbulent  flow  with  normal 
shock  stagnation  pressure  and  with  two  pressure  distributions,  corres¬ 
ponding  to  N  £  1/4  and  for  N  =  1.  13.  The  same  method  was  also 
employed  fo£  the  pressure  distribution  corresponding  to  N  <  1/3  but 
with  oblique  shock  stagnation  pressure.  3  ~ 

• 

Consideration  of  Figure  10  indicates  that  most  of  the  data  obtained 

with  nitrogen  involve  transitional  and  turbulent  flow.  The  heating  rates 

over  the  front  portion  of  the  cone  are  reduced  by  the  mass  transfer  but  for 

x;>  5  the  effect  of  mass  transfer  is  to  cause  earlier  transition  and  therefore 

heating  rates  whj.ch  are  greater  with  than  without  mass  transfer.  The 

maximum  heating  rate's  obtained  experimentally  are  greater  than  those 

predicted  by  the  F.  P.  R.  E.  method  when  the  stagnation  pressure  of  the  flow 

external  to  the  boundary  layer  corresponds  to  that  behind  a  normal  shock. 

However,  if  the  stagnation  pressure  is  assumed  to  be  that  behind  the 

oblique  shock,  the  predicted  turbulent  heating  without  mass  transfer  is 

greater  than  determined  experimentally.  The  results  in  reference  15  for 

the  interaction  of  a.  turbulent  boundary  layer*with  the  vortical  layer  induced 

by  the  bow  shock  indicate  that  some  influence  may  be  expected  for  x  ;>  10 

under  the  test  conditions  of  ‘Figure  10  with  no  mass  transfer.  Therefore, 

the  observed  rates  with  mass  transfer  may  involve  vortical  interaction. 

In  view  of  this  complication,* it  is  not  possible  to  ascertain  the  influence 

of  the  mass* transfer  on  the  heat  transfer  after  transition-has  occurred. 

•  • 

Transition  Behavior 


From  a  representation  of  the  heat  transfer  data  for  nitrogen 
in  the  form  shown  in  Figure  10  it  is  possible  to  estimate  the  location 
of  transition  for  various  rates  of  injection  and  various  Reynolds 
numbers.  The  results  are  shown  in  Figure  11  for  the  test  data  corres¬ 
ponding  to  N  £  1/3;  as  seen  qualitatively  from  Figure  10  and  discussed 
previously,  injection  results  in  transition  at  a  smaller  value  of  x.  It  can 
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also  be  seen  from  Figure  11  that  transition  moves  downstream  (x  increasing) 
as  the  Reynolds  number  decreases. 

It  is  customary  to  consider  transitioh  behavior  in  terms  of  a  Reynolds 
•  number  based  on  momentum  thickness  (cf.  for  example  reference  14).. 
However,  the  analysis  of  the  data  of  Figure  11  in  these  terms  awaits  a  theor¬ 
etical  treatment  of  the  downstream  behavior  of  the  laminar  boundary  layer. 
Presently,  the  data  of  Figure  11  must  be  considered  only  indicative.  ' 


IV.  CONCLUSIONS 


•  There  has  been  carried  out  a.n  investigation  of  the  influence  of 
.localized  mass  transfer  at  the  nose  of  a  slender  cone.  Two  gaseous 
coolants  were  considered,  nitrogen  and  helium.  The  experimental 
research  involved  tests  in  the  Mach  number  8.  0  wind  tunnel  at  PIBAL. 

The  rate  of  coolant  mass  flow  and  the  Reynolds  number  of  the  main  stream 
-were  varied  widely.  Measurements  were  made  of  the  pressure  and  heat 
transfer  distributions;  schlieren  photographs  were  taken. 

Three  effects  of  the  mass  transfer  are  examined:  the  alteration  of 
shock  shape  and  pressure  distribution,  of  heat  transfer  to  the  downstream 
conical  surface,  and  of  transition.  With  respect  to  the  first  effect  a 
the'oretical  analysis  of  the  influence  of  mass  transfer  on  the  shock  stand¬ 
off  distance  has  been  carried  out.  The  model  for  the  shock  layer  proposed 
by  Lighthill  is  combined  with  a  model  suggested  by  the  results  of  reference 
2  for  the  viscous  layer  formed  by  the  coolant.  Accordingly, the  flow  at  the 

•  nose  is  characterized  by  two  incompressible  rotational  layers  with  a 
discontinuity  in  tangential  velocity  at  the  interface.  This  analysis  which  is 
applicable  to  both  nitrogen  and  helium  suggests  an  inviscid  flow  parameter 
N  which  is  related  to  the  mass  transfer  and  which  is  shown  to  correlate 

_  o  • 

.  the  observed  alteration  in  shock  stand-off  distance  by  coolant  injection. 

The  mass  transfer  is  found  to  have  a  negligible  effect  on  the 
pressure  distribution. on  the  cone  for  combinations  of  mass  transfer  and 
Reynolds  numbers  resulting  in  N  <1/3.  Thus  for  flows  satisfying  this 
inequality  laminar  boundary  laye?-  effects  are  correlated  by  the  parameters 
presented  in  reference  1 ,  i.  e.  ,  Nx  and  N^,  For  flows  corresponding  to 

N  >  1/3  the  pressure  distribution  is  altered  by  the  mass  transfer  so  that 
no  extrapolation  of  results  to  flow  conditions  other  than  those  of  the  tests 
is  possible. 

•  *  •  •  • 

With  respect  to  the  influence  of  mass  transfer  on  heat  transfer  it 
is  found  as  shown  previously  that  heliiim  is  the  superior  coolant.  With 
helium  at  relatively  low  rates  of  coolant  flow  the  peak  heating  occurs . 
downstream  on  the  cone  and  is  an  order  of  magnitude  less  than  would 
occur  at  the  stagnation  point  if  there  were  no  mass  transfer. 


35 


The  injection  of  nitrogen  is  found  to  lead  to  earlier  transition.  Thus 
on  the  cone  the  peak  heating  with  injection  is  greater  than  without  injection 
although  considerably  less  than  at  the  stagnation  point.  The  behavior  of 
transition  with  injection  is  deduced  from  the  heat  transfer  results;  as  might 
be  expected  transition  moves  forward  on  the  cone  with  increasing  mass 
transfer  and  increasing  Reynolds  number. 
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APPENDIX  I  -  THE  SHOCK  STAND-OFF  DISTANCE 
WITH  STAGNATION  POINT  MASS  TRANSFER 


The  shock  stand-off  distance  has  been  determined  for  flow  about 
a  sphere  by  Lighthill  (reference  16)  and  about  a  circular  cylinder  by 
Whitham  (reference  17).  Both  analyses  utilize  the  same  approximations. 
These  consist  of: 

( 1)  applying  strong  shock  conditions  which  are  valid  for  large 
free- stream  Mach  number; 

(2)  considering  the  flow  field  between  the  body  and  the  shock  to 
be  incompressible  and  rotational;  and 

(3)  assuming  the  shock  to  be  a  sphere  or  a  circular  cylinder.' 

It  is  generally  accepted  that  these  solutions  accux-ately  predict  the  shock 
stand-off  distance  but  not  the  pressure  distribution. 

The  purpose  of  the  present  analysis  is  to  extend  the  results  of 
Lighthill  to  include  the  effect  of  mass  transfer  from  a  spherical  body.  Note 
that  the  injection  region  is  assumed  to  extend  around  the  entire  forward 
portion  of  the  body.  The  treatment  of  a  c-ircular  cylinder  as  presented 
by  Whitham  can  also  be  extended  to  include  mass  transfer  effects. 

The  flow  is  shown  schematically  in  Figure  12.  A  spherical  shock 

of  radius  r  =  A  +  R  is  associated  with  a  sphere  of  radius  R  .  The  flow 
so  r  o 

between  the  shock  and  the  body  is  divided  into  two  regions  by  an  interface 
of  radius  r  .  In  the  outer  layer  between  the  shock  and  the  interface  the 
flow  is  of  constant  density  and  is  rotational  because  of  the  shock  curvature. 
The  inner  layer  is  between  the  interface  and  the  body;  the  gas  therein  is 
also  incompressible  but,  in  general,  of  compositiorf  and  temperature  which 
are  fixed  but  different  from  that  in  the  outer  layer.  The  flow  in  the  inner 
layer  *is  also  rotational.  *  At  the  interface  the  radial  velocity  is  zero; 
across  the  interface  the  pressure  is  continuous  but  the  vorticity  and  tan¬ 
gential  velocity  are  discontinuous. 

The  equations  of  motion  for  an  axisymmetric ,  incompressible, 
rotational  flow  in  spherical  coordinates  reduce  to: 

—  [rvuj  +  ~[vu]=0  ( I-  L) 

0r‘  r  8cp  ® 


The  rotationality  in  the  inner  layer  must  be  considered  to  be  due  to  the 
shear  stresses  which  actually  occur  between  the  inner  and  outer  flows  and 
which  disappear  in  the  idealization  of  the  flow  into  two  inviscid  but 
rotational  regions, 
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Introduce  the  Stokes  stream,  function,  which  satisfies  mass  con¬ 
tinuity,  and  which  is  related  to  the  radial  and  tangential  velocity  compon¬ 
ents  by: 


v  = 
cp 


r  sin  cp  3r 


v  =  - 
r 


1 


ra  sin  cp  8cp 

The  vorticity  in  terms  of  the  stream  function  is  given  by: 


to  = 


1 '  ra  ®!l  +  ®!A  -  1*  cot  cp  ] 
aVa  acra  dcp 


r  sin  cp 


(1-2) 

(1-3) 


d-4) 


Outer  Layer 

Lighthill's  solution  applies  for  the  layer  between  the  spherical 

shock  (r=r  )  and  the  interface  (r=r  ).  Accordingly,  the  velocity  com- 
s  o 

ponents  at  the  outer  surface  of  the  interface  are  given  by: 


v 

r 


=  0 


v 

Cp 


=  (v  ) 
cp  e 


V  sin  cp 

— - r(K-l)3  (r  J"a  -K(K-4)i 

2K  8 


(1-5) 


where  K  is  the  static  density  ratio  across  a  normal  shock;  i.  e.  ,  Pg/p^* 

r  =  r/r  ,  and  V  is  the  free  stream  velocity.  The  location  of  the  shock 
o  oo 

relative  to  the  interface  is  determined  by  solving  the  following  fifth  degree 
algebraic  equation  for  r  =  r5/r0- 


2(K-6)(?  )B  -  IMjSlll  r  3  +  3(K-  1)  =  0 

®  XT  1  8 


(1-6) 


If  there  is  no  inner  layer,  i.  e,  ,  no  injection  this  equation  with  rQ-*  Rq 

gives  the  shock  stand-off  distance  A/r  =  a/R  =  r  -1. 
e  oos 

Conditions  at  the  interface 


Now,  applying  the  incompressible  Bernoulli  equation  on  both  sides 
of  the  interface,  and  equating  the  static  pressures  all  along  the  interface 
yields:  - 


41 


42 


where* 


V  f  0  W  i 

P  =^.J(K-l)a  (r  )“»  -  K(K-4)^  l]* 


(1-14) 


and  where  R  s  R  /r  . 

o  o 

The  radial  velocity  at  the  body  surface  can  be  determined  from  the 
stream  function  and  is  found  to  be: 

2p  f  (2  +  R  "a)(2R"a  -  5  +  3R  af| 

(v  )  =  — cos  cpi(R  )'3-  1+ - *. - °  ■  (I-  15) 

r  v  3  12  R  a  -  10  -  2R  _a 

L  o  o  J 


It  is  noted  that  the  solution  corresponds  to  a  non-uniform  injection 
rate  over  the  surface  according  to  cos  cp-l-cpP/2  +  .... 

The  mass  transfer  rate  of  interest  for  comparison  with  experiment 
is  obtained  by  integrating  the  injection  velocity  over  the  porous  region 
and  becomes : 


T  W  ll 

_  o  C  3  . 

m  =  irR  q  - - -  si 

c  °  ‘  e  0  W 

I  w  e  J 


inacp.  — —  [ (K  -  l)a  x  ~a  -  KtK-4)  ] 
1  3K  8 


(R  ~s  +  2)(2R  ~a  -  5  +  3R  a) 


Utilizing  strong  shock  conditions,  the  free  stream  velocity  (V  )  can  be 
related  to  stagnation  conditions  behind  the  normal  shock  and  tSe  density 
ratio  (K)  by. 


It  should  be  noted  that  p  is  the  stagnation  point  velocity  gradient 
dvcp/dcp]c  at  the  interfaci.  IE  the  velocity  gradient  of  the  external  flow 

(at  the  interface  is  denoted  by  p  ,  the  relation  between  the  two  gradients 

®  0  W  | 

r  y?  Q  3 

across  the  interface  is  given  by:  p  =  p  [ - J  . 


There  results  the  non-dimensional  similarity  parameter,  N  given  by: 


m 


N  = 


R  3(p  pa  )3 
o  rs  Ks 
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-  sin8  cp.  [(K-  l)a  r/3  -  K(K~4)  ] 
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R  s  -  1  +  - 2 - 2__ 
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12R  3  -  10  -  2R  "3 
o  o 


(1-18) 


Consider  the  application  of  this  analysis  to  a  particular  example. 

If  M  and  the  other  free  stream  quantities  are  known,  K  is  determined  by 
normal  shock  conditions.  The  parameter  N  and  the  angle  cp.  are  assumed 

known.  Noyr  Eq.  (1-6)  can  be  solved  for  J  and  Eq._(I-l8)  can  then  be 

solved  for  R  .  Note  that  A/R  =  (r  /R  )-Sl=(r/R  -1.  In  Figure  6  the 

o  o  &  o  so 

variation  of  a/R  with  N  for  M  =8^0  has  been  given  and  compared  to 
.  .  o  a  oo 

experiment. 


Since  the  stagnation  enthalpy  achieved  in  the  wind  tunnel  is  much 
less  than  that  corresponding  to  free  flight  conditions  at  M  =8.  0,  the  strong 
shock  approximations  are  not  strictly  valid  for  the  analysis  of  the  experi¬ 
mental  data.  As  a  result,  in  computing  the  theoretical  curve  presented  in 
Figure  6,  the  approximation  attendant  with  Eq.  (1-17)  was  not  used  so  that 
Eq.  ( X-  18)  was  modified  somewhat. 


Tangential  velocity  distribution 


It  is  further  of  interest  to  consider  the  variation  of  the  tangential 
velocity  distribution  with  N  and  shown  in  Figure  13.  It  is  noted  that  for 
relatively  small  values  of  3N  the  ratio  (r  /R  )  ~  1,  and  the  tangential 
velocity  distribution  is  close  3to  linear;  thus  tSe  inviscid  model  for  the 
stagnation  point  boundary  layer  gives  in  this  case  the  same  tangential 
velocity  distribution  as  the  boundary  layer  calculations  which  were  given 
in  reference  2  and  which  were  based  on  6/R  <  <  1.  However,  for  increasing 
N  ,  the  ratio  r  /R  increases  and  sphericity  results  in  a  nonlinear  dis¬ 
tribution  of  tangential  velocity.  Presumably,  the  boundary  layer  equations 
with  6/R  ~  1  would  result  in  a  similar  distribution  except  in  a  relatively 

thin  regi8n  which  is  close  to  the  outer  edge  and  which  involves  significant 
shear  stresses. 

In  the  outer  flow  layer,  with  the  assumption  of  a  spherical  shock, 
the  ratio  of  vorticity  to  distance  from  the  center  line  is  constant  throughout 
the  region.  For  the  boundary  conditions  utilized  in  the  solution  of  the  inner 
flow  layer,  it  is  noted  that  this  ratio  is  also  constant  in  this  flow  region. 
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FIG,  (13)  TANGENTIAL  VELOCITY  PROFILES  FOR  VARIOUS 
MASS  TRANSFER  RATES  (NITROGEN) 


Empirical  correction  for  velocity  gradient 

B:  is  well  known  that  the  velocity  gradient  §  is  not  well  predicted 
by  the  Xiigb.th.ill  theory,  the  error  being  due  to  the  f  seumed  spherical 
shockwave.  An  empirical  correction  to  the  analysis  presented  her®  can, 
however,  be  madej  if  8  is  computed  from  Newtonian  theory  applied  at 
the  interface,  then  Eq,  (1-18)  becomes 


N 
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2irv/  2 
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sina  cp. 
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(2  4  R  “a)(2R  _a  -  5'4  3R  sji 

1  ^ - 2 - - 2 - - - i_A 

I2R  3  -  10  -  2R  :a  .  .  • 

O  O  J  : . 

(I- 18a) 


For  M^sS.  0  the  results  of  the  determination  of  A/R  from 
Eq.  (I- 18a)  are  also  shown  in  Figure  6  and  are  seen  to  be  in  somewhat 
better  agreement  with  experiment. 
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